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In  this  study,  pyrolysis  kinetics  of  periphytic  microalgae  consortia  grown  using  swine  manure  slurry  in 
two  seasonal  climatic  patterns  in  northwest  Arkansas  were  investigated.  Four  heating  rates  (5,  10,  20 
and  40  °C  min  1 )  were  used  to  determine  the  pyrolysis  kinetics.  Differences  in  proximate,  ultimate, 
and  heating  value  analyses  reflected  variability  in  growing  substrate  conditions,  i.e.,  flocculant  use, 
manure  slurry  dilution,  and  differences  in  diurnal  solar  radiation  and  air  temperature  regimes.  Peak 
decomposition  temperature  in  algal  harvests  varied  with  changing  the  heating  rate.  Analyzing  pyrolysis 
kinetics  using  differential  and  integral  isoconversional  methods  (Friedman,  Flynn-Wall-Ozawa,  and 
Kissinger-Akahira-Sunose)  showed  strong  dependency  of  apparent  activation  energy  on  the  degree  of 
conversion  suggesting  parallel  reaction  scheme.  Consequently,  the  weight  loss  data  in  each  thermogravi- 
metric  test  was  modeled  using  independent  parallel  reactions  (IPR).  The  quality  of  fit  (QOF)  for  the  model 
ranged  between  2.09%  and  3.31%  indicating  a  good  agreement  with  the  experimental  data. 

©  2014  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

Thermogravimetry  is  a  powerful  and  versatile  tool  in  under¬ 
standing  and  modeling  biomass  reactions.  It  is  also  a  quick  and 
reliable  approach  to  determine  the  moisture,  organic  matter,  and 
ash  contents  in  biomass  (Garcia  et  al„  2013).  Other  uses  of  thermo- 
gravimetric  analysis  (TGA)  include  quantification  of  the  hemicellu- 
lose  and  a-cellulose  contents  in  wood  (Carrier  et  al.,  2011),  as  well 
as  evaluation  of  digestate  stability  during  anaerobic  digestion 
(Gomez  et  al.,  2007).  The  predominant  use  of  TGA,  however,  is  to 
determine  the  rate  of  thermal  decomposition  for  various  feed¬ 
stocks.  A  survey  of  the  literature  shows  the  importance  of  this 
analysis  to  the  determination  of  decomposition  kinetics  of  woody 
biomass  (Gomez  et  al.,  2005),  crop  residue  (Mani  et  al.,  2010), 
animal  manures  (Sharara  and  Sadaka,  2014),  and  municipal  solid 
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waste  (Sanchez  et  al.,  2009).  However,  the  recent  interest  in 
aquatic  biomass,  as  a  biofuel  and  bioenergy  crop,  brought  to  light 
the  dearth  of  thermogravimetric  studies  on  this  class  of  biomass. 

Historically,  microalgae  mass-production  was  mostly  accom¬ 
plished  to  support  aquaculture  systems  (Muller-Feuga,  2000)  as 
well  as  to  extract  important  bioactive  compounds  (Skulberg, 
2004).  Recently,  however,  interest  in  sustainable  biofuel  sources 
brought  attention  to  microalgae  due  to  the  high  lipid  content  in 
certain  species,  which  can  be  converted  via  “transesterification” 
to  biodiesel.  The  oil  content  in  microalgae  species  such  as  Schizo- 
chytrium  sp.  can  exceed  75%  on  a  dry-weight  basis  (Demirbas 
and  Fatih  Demirbas,  201 1 ).  A  vital  ecological  service  that  microal¬ 
gae  could  provide  is  treatment  of  nutrient-rich  wastewater  efflu¬ 
ents,  i.e.,  phycoremediation.  This  particular  service  is  becoming 
increasingly  crucial  with  the  increase  in  global  urbanization,  indus¬ 
trialization,  and  intensive  cropping  activities  (Heathwaite,  2010). 

Contamination  of  surface  and  ground  water  with  nitrogen  (N), 
phosphorous  (P),  polycyclic  aromatic  hydrocarbons  (PAH),  and 
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heavy  metals  is  the  primary  cause  behind  hypoxic  conditions, 
eutrophication,  and  the  collapse  of  aquatic  habitats  in  rivers  and 
lakes  (Conley  et  al.,  2009).  Uptake  and  biosorption  of  these  con¬ 
taminants  by  algal  species  has  been  investigated  on  a  variety  of 
wastewater  substrates,  i.e„  municipal,  industrial,  and  agricultural 
wastewater  streams  (Pittman  et  al.,  201 1 ).  Adopting  phycoremedi- 
ation  in  livestock  wastewater  management  can  offer  an 
added-advantage  by  minimizing  phosphorus  loss.  Algal  species 
were  successfully  grown  on  raw  swine  manure  effluents, 
9  g  m-2  day-1,  but  were  found  to  be  sensitive  to  the  loading  rates 
of  N  and  P  (Kebede-Westhead  et  al.,  2006). 

The  harvested  phycoremediation  algae  can  be  directed  to  a  vari¬ 
ety  of  uses  such  as  composting,  anaerobic  digestion,  and  the  extrac¬ 
tion  of  lipids  and  sugars  for  biodiesel  and  ethanol  production. 
Alternatively,  the  algal  biomass  can  be  utilized  through  thermo¬ 
chemical  conversion  processes.  The  end-use  for  a  specific  algae  har¬ 
vest  is  highly  dependent  on  its  composition.  Algae  used  in  water 
remediation  are  typically  indigenous  species  grown  in  open  systems 
that  results  in  a  diverse  species  consortium  in  the  collected  biomass. 
Thermogravimetric  analyses  can  assist  in  the  effort  to  characterize 
the  harvested  microalgae  and  to  direct  them  to  the  optimal  end- 
uses.  For  instance,  TGA  decomposition  rate  curves  were  proven 
effective  in  detecting  and  quantifying  the  lipid  contents  in  Chlorella 
sp.  (Na  et  al.,  2011).  Understanding  the  decomposition  kinetics  of 
microalgae  in  various  atmospheres,  in  addition  to  providing  insight 
into  the  proximate  composition,  is  essential  to  the  design  and  opti¬ 
mization  of  thermochemical  conversion  processes.  There  is  a  need  in 
research  for  studies  covering  the  thermal  decomposition  kinetics  of 
mixed  algae  consortia  grown  for  water  remediation  purposes,  espe¬ 
cially  as  part  of  sustainable  livestock  production  scenarios.  The  goal 
of  this  research  is  to  study  the  thermal  decomposition  of  indigenous 
periphytic  microalgae  grown  in  an  open  water-remediation  system 
for  a  swine  production  facility  under  nitrogen  environment  using 
thermogravimetric  analysis  (TGA). 


2.  Methods 

2.3.  Swine  wastewater  treatment  system 

Algal  biomass  was  produced  using  open  channel  raceways  lined 
with  a  growth  medium  for  periphyton  attachment  and  irrigated  by 
circulation  of  swine-manure  based  wastewater.  Algae  consortia 
produced  using  two  systems  were  investigated  in  this  study.  Sys¬ 
tem  1  was  a  small  system  with  8  raceways  that  were  each  15  cm 
wide  by  3.0  m  long.  It  was  sited  outdoors  adjacent  to  the  University 
of  Arkansas  Biological  and  Agricultural  Engineering  Research  Lab, 
Fayetteville,  Arkansas.  System  2  was  a  larger  pilot  system  with  4 
raceways  that  were  each  1.5  m  wide  by  61  m  long.  It  was  con¬ 
structed  at  the  University  of  Arkansas  Swine  Research  Grower 
Facility,  near  Savoy  Arkansas.  The  raceways  in  both  systems  were 
lined  with  specialty  fibrous  carpeting  (proprietary  fiber  selection 
and  layout  intended  to  optimize  periphyton  growth,  provided  by 
Interface,  Inc.)  to  maximize  algae  attachment  during  water  circula¬ 
tion.  System  1  was  operated  during  the  summer  of  2013  (June- 
July)  on  swine  manure  slurries  at  different  degrees  of  dilution  to 
vary  targeted  starting  ammonia  (NH3)  concentrations  between  5 
and  40  mg  L-1  (mg  =  10-6  kg).  Details  of  the  NH3  loading  influence 
will  be  discussed  in  a  separate  publication.  System  1  was  seeded 
with  mixed  consortia  which  were  collected  from  a  local  stream 
in  Fayetteville,  AR  whereas  System  2  was  seeded  from  the  mixed 
consortia  harvested  from  System  1.  Wastewater  used  in  System  2 
was  pumped  from  swine  slurry  storage  lagoon,  then  treated  with 
alum  (aluminum  sulfate  14-hydrate  granules,  A12(S04)314H20, 
addition  rate  2gL-1)  to  flocculate  solids.  After  a  24  h  settling 
period,  the  undiluted  supernatant  was  removed  for  circulation  in 


the  raceway.  System  2  was  operated  during  the  fall  of  2013 
(November-December).  Both  systems  were  open  during  produc¬ 
tion  without  any  control  over  the  growing  algal  species. 

2.2.  Algae  collection  and  preparation  methods 

Harvest  of  algae  was  accomplished  manually,  in  5-day  cycles, 
using  a  rubber-bladed  squeegee  for  removal  and  collection  of 
attached  algae.  Harvested  algae  were  dried  immediately  or  else 
stored  at  4  °C  for  24  h  before  drying.  Subsamples  of  fresh  algae 
were  stored  for  analysis  and  species  identification.  One  composite 
sample  was  collected  from  each  algae  consortium  for  the  thermo¬ 
gravimetric  analysis  (TGA).  Composite  algae  samples  were  dried  at 
105  °C  for  24  h  before  grinding  to  1  mm  (Thomas  Wiley®  cutting 
mill-Model  3383L10,  Swedesboro,  NJ).  A  second  grinding  step, 
using  a  cutting  mill  (Polymix  PX-MFC  90  D,  Kinematica  AG, 
Switzerland),  was  added  to  homogenize  the  dried  samples  and 
minimize  mass  transport  resistance  during  thermogravimetric 
analysis.  The  algal  solids  used  in  the  following  tests  and  analyses 
all  passed  through  a  200  pm  sieve  (pm  =  10-6  m). 

2.3.  Proximate,  ultimate,  heating  value  and  pH  analyses 

The  remaining  moisture  was  determined,  in  triplet  samples,  as 
the  weight  loss  after  drying  at  105  °C,  ASTM  E871  -  82  (2006). 
Standard  methods  were  also  used  to  determine  the  volatile  matter, 
ASTM  E872  -  82  (2006),  and  ash  content,  ASTM  D2974  -  87  (2007), 
while  the  fixed  carbon  (%)  was  determined  by  difference.  Complete 
elemental  analyses  of  representative  samples  from  each  algae 
harvest  were  performed  in  a  specialized  diagnostic  laboratory 
(Huffman  laboratories,  Golden,  CO,  USA).  The  heating  values  were 
determined  on  sample  triplets,  according  to  standard  ASTM  D5865 
-  12  (2012),  using  bomb  calorimetry  (Parr  instruments,  Model 
1341,  Moline,  IL,  USA).  The  pH  of  algal  biomass  was  determined 
using  a  pH  electrode  (SB70P,  SympHony,  VWR,  Radnor,  PA,  USA) 
after  the  dry,  ground  samples  were  diluted  with  deionized  water, 
10  mL  per  1  g  of  sample,  then  vigorously  stirred  and  allowed  to 
stand  for  1  h  before  measurement. 

2.4.  Thermogravimetric  analysis  methodology 

A  thermogravimetric  analyzer  (Model  TGA  4000,  PerkinElmer, 
Inc.  Waltham,  MA,  USA)  was  used  to  study  the  decomposition 
behavior  of  the  two  algal  harvests.  Prior  to  the  algae  decomposition 
tests,  a  curie-point  temperature  calibration  was  performed  using 
three  reference  materials,  i.e.,  alumel,  perkalloy,  and  Iron,  accord¬ 
ing  to  the  manufacturer’s  guidelines.  Similarly,  a  weight  calibration 
was  performed,  before  tests,  using  the  manufacturer  supplied 
reference  weight.  Pyrolysis  of  each  of  two  algae  biomasses  was 
studied  under  four  different  heating  rates  (5,  10,  20,  and 
40°Cmin-1)  from  30  to  800  °C.  Nitrogen  gas  was  used  to  purge 
the  sample  (30mLmin-1)  to  simulate  pyrolysis  conditions.  The 
sample  size  was  consistently  kept  at  5  ±  0.5  mg  to  minimize  devi¬ 
ation  (lag)  between  measured  and  actual  sample  temperatures, 
and  to  also  ensure  that  the  decomposition  reactions  were  not 
transport-limited.  For  each  sample,  blank,  clean  crucibles  were 
tested  using  the  same  thermal  decomposition  programs  in  order 
to  adjust  the  weight  baseline  by  compensating  for  the  drag  force 
acting  on  the  crucible. 

2.5.  Decomposition  kinetics 

The  decomposition  is  often  expressed  in  terms  of  the  conver¬ 
sion  (a)  which  describes  the  change  in  sample  weight,  in  a  dimen¬ 
sionless  form. 
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W0,  Wt,  and  W^  are  the  sample  weights  at  the  beginning,  at  time  t, 
and  at  the  end  of  the  decomposition  stage,  respectively.  The  rate  of 
conversion  (da/dt)  is  often  expressed  using  an  Arrhenius  type 
expression. 

W  =  AeXp(w)f(a)  (2) 


2.5. 3.2.  Kissinger  method.  This  method  determines  the  activation 
energy  (Ea)  and  the  pre-exponential  factor  (A)  using  the  tempera¬ 
tures,  Tmax  values,  that  correspond  to  peak  decomposition  rates, 
(da/dt)max  in  thermogravimetric  tests  under  different  heating  rates 
(p).  The  temperature  integral  in  Eq.  (4)  does  not  yield  an  analytical, 
closed-form  solution.  Therefore,  it  is  alternatively  presented  as 
follows  (Starink,  2003) 

AF 

g(a)=1£p(x)  (6) 


A,  £a,  R,  and  T  are  the  pre-exponential  coefficient  (frequency  factor), 
the  activation  energy,  the  universal  gas  constant,  and  the  sample 
absolute  temperature,  respectively.  f[ a)  represents  the  kinetic 
model  that  describes  the  rate  of  conversion  dependence  on  the 
conversion,  e.g.,  an  n-order  reaction  model:  /(a)  =  (1  -  a)". 

Under  a  constant  heating-rate  {ft  =  dT/dt),  the  time-dependence 
of  the  conversion  rate  can  be  transformed  to  a  temperature- 
dependence  which  can  be  used  to  rewrite  the  differential  form 
(Eq.  (3))  or  the  integral  form  (Eq.  (4))  of  the  decomposition  kinetic 
expression. 

S=^exp©/(a)  (3) 


where  x  =  E^RT,  and  the  function  p{x)  represents 

/*o°  e-x 

p(x)  =  J  Y2 ~dx  (7) 

Many  approximations  were  derived  to  determine  p(x )  numeri¬ 
cally.  In  the  Kissinger  method,  the  Murray  and  White  approxima¬ 
tion  (Starink,  2003),  p(x)  =  e~*/x2,  is  used.  The  underlying 
assumption  in  this  method  is  that  the  decomposition  follows  a 
first-order  reaction  model,  i.e.,  /(a)  =  (l  -  a).  Eq.  (8)  details  the 
relationship  between  Tmax,  p,  and  £a  and  A. 


g(a)=r?exp©dr  w 

The  following  sections  will  detail  the  different  model-free 
(isoconversional)  methods  used  to  determine  the  apparent  activa¬ 
tion  energies  for  algae  pyrolysis. 

2.5.3.  Model-free  (isoconversional)  methods 

Isoconversional  methods  overcome  the  requirement  of  deter¬ 
mining  the  reaction  model, /(a),  in  order  to  determine  the  activa¬ 
tion  energy,  £a.  This  is  accomplished  by  simultaneously  analyzing 
decomposition  curves  generated  under  different  heating  rates  to 
extract  the  apparent  kinetics  data,  i.e.,  £a  and  In  A,  corresponding 
to  each  degree  of  conversion  (a).  Isoconversional  methods  are  pop¬ 
ular  in  biomass  decomposition  studies  due  to  the  fact  that  biomass, 
a  natural  biopolymer,  undergoes  a  series  of  overlapping  reactions, 
during  pyrolysis  or  oxidation,  which  cannot  be  described  accu¬ 
rately  by  one-step,  global  reaction  model.  The  downside  in  isocon¬ 
versional  methods,  however,  is  the  inability  to  straightforwardly 
determine  an  exact  model  expression, /(a)  or  g(ot),  to  describe  the 
entire  decomposition. 

In  this  research,  one  differential  isoconversional  method, 
Friedman  method  (Friedman,  1965),  and  three  integral  isoconver¬ 
sional  methods:  Kissinger’s  (Blaine  and  Kissinger,  2012), 
Flynn-Wall-Ozawa’s  (Flynn  and  Wall,  1966;  Ozawa,  1965),  and 

Kissinger-Akahira-Sunose’s  (Starink,  2003)  were  used  to  deter¬ 
mine  the  pyrolysis  kinetics  of  the  two  algae  harvests  under  study. 
Below  is  a  brief  description  of  each  the  methods. 

2.5. 3. 3.  Friedman  method.  Re-arranging  Eq.  (3),  then  linearizing  it 
by  taking  the  natural  log  of  both  sides  of  the  equation 

In  =ln(A/(a))-|i  (5) 


By  plotting  In  P/T^x  against  (1/Tmax),  the  activation  energy 
(£a/R)  and  the  pre-exponential  factor,  In  (AR/£a),  terms  can  be 
determined,  respectively,  as  the  slope  and  the  intercept  of  the 
resulting  straight-line. 


2. 5. 1.3.  Kissinger-Akahira-Sunose  (I<AS)  method.  A  modified  form  of 
Kissinger’s  method  that  was  described  earlier,  KAS  method  dis¬ 
penses  with  the  kinetic  model  assumption,  and  the  use  of  only 
the  peak  decomposition  data.  Instead,  KAS  is  an  integral  method 
that  relies  on  the  different  (a)  values  instead  of  the  single  value 
corresponding  to  peak  decomposition  (da/dt). 

Eq.  (8)  can  then  be  rearranged  and  linearized 


In  (g(oe)) 


Ea 

RT 


(9) 


The  slope  of  the  straight  lines  formed  by  plotting  In  (PIT2) 
against  1/T  at  each  degree  of  conversion  (a),  with  the  points 
corresponding  to  the  different  p  values,  yields  the  activation 
energy,  £a,  corresponding  to  that  conversion  degree. 


2.5. 3.4.  Flynn-Wall-Ozawa  (FWO)  method.  In  this  method,  the 
Doyle’s  approximation  (Doyle,  1962)  to  the  temperature  integral, 
p(x),  is  applied:  p(x)  s  -5.33  -  1.05x.  The  linearized  form  of 
Eq.  (4)  becomes: 


In  p  =  In 


A£a 

Rg(a) 


-  5.33  -  1.052 


(10) 


At  each  degree  of  conversion  (a),  In  p  and  1/T  corresponding  to 
each  heating  rate  are  fitted  into  a  straight  line.  The  slope, 
(-1.052 £a/R),  represents  the  apparent  activation  energy  term, 
while  the  intercept  is  a  coupled  expression  of  the  reaction  model 
in  the  integral  form,  g(oc),  the  apparent  activation  energy  £a,  and 
the  frequency  factor,  A. 


At  each  oe,  the  above  equation  describes  the  linear  relationship 
between  1/T  and  In  (p  doc/dT)  with  each  point  representing  a  tested 
heating  rate.  The  slope  of  this  line  represents  the  activation  energy 
term,  EJRT,  at  this  conversion  degree  whereas  the  intercept 
represents  an  expression  combining  the  reaction  model  and  the 
pre-exponential  factor,  In  (Af{  a)).  Calculating  the  slopes  and  inter¬ 
cepts  at  different  a  values,  between  0.05  and  0.90,  describe  the 
kinetics  of  decomposition  as  a  function  of  the  conversion  degree. 


2.5.2.  Kinetic  modeling  of  algae  pyrolysis 

In  order  to  investigate  the  pyrolysis  of  the  two  algae  consortia, 
the  weight  loss  was  modeled  as  a  series  of  overlapping  indepen¬ 
dent,  parallel  reactions  taking  place  between  100  and  700  °C.  This 
approach  was  presented  elsewhere  as  a  deconvolution  step  neces¬ 
sary  to  understand  the  devolatilization  of  a  biopolymer  by  model¬ 
ing  the  various  pseudo-components  contributing  to  its  overall 
decomposition.  In  this  study,  algae  pyrolysis  was  modeled  as  4 
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independent,  1st  order  reactions.  The  rate  of  conversion  (da/dt)  for 
pseudo-component  (N)  is: 

^  =  ANexp(-^j(\-<xN),  whereN  =  1 . .  .4  (11) 

Consequently,  the  total  rate  of  conversion  ( dxTotai/dt )  can  be 
determined  from  the  individual  pseudo-components  using  the 
following  expression: 


t^Total 

dt 


EW- 


dCLft 

dt 


N 

andJ^W;  =  1 

i-l 


(12) 


The  quantity,  w„  represents  the  contribution  of  each  individual 
pseudo-component  to  the  overall  sample  conversion.  The  conver¬ 
sion  rate  can  be  converted  to  the  mass  loss  rate  using  the  following 
relation: 


dmcalc  /m  rr,  \  ^aT0tal 

dt  dt 


(13) 


mQ  and  mf  represent  the  initial  and  final  sample  weights.  Numerical 
integration  between  the  temperatures  of  interest,  i.e.,  100-700  °C  in 
this  study,  yields  the  computed  overall  weight,  mcalc.  The  parame¬ 
ters  of  each  reaction,  i.e.,  A,  E3,  and  w  were  determined  by  nonlinear 
minimization  of  an  objective  function  (O.F.)  which  represents  the 
sum  of  squared  differences  between  observed  sample  mass  loss 
rate,  dinobsldt,  and  the  calculated  mass  loss  rate,  dmc3ijdt: 


(14) 


In  this  study,  the  total  number  of  points,  Np,  used  in  reaction 
modeling  for  each  thermogram  was  300  points.  The  quality  of  fit, 
QOF  (%),  for  the  weight  loss  derivative  (DTG)  and  the  weight  loss 
(TG)  were  determined  using  the  following  expressions: 


QOFdtg(%)  =  100 


I  ^ 

l  * 

dt  ) 

/ 

dmobs 

V 

Np 

/ 

dt 

(15) 


3.2.  Proximate  analysis 

Table  1  shows  the  results  of  proximate  analysis  of  the  two  algae. 
The  low  moisture  content  in  both  harvests,  around  5%  wet-basis, 
represents  the  moisture  absorbed  after  drying  and  during  sample 
preparation.  Comparison  of  proximate  analysis  results  for  the 
two  algae  harvests  showed  significant  differences  (t-test, 
p<0.01)  in  volatile  matter,  fixed  carbon,  and  ash  contents.  The 
ash  contents  observed  in  both  algae  consortia:  20%  in  Algae  1 
and  32%  in  Algae  2  are  noticeably  high  in  comparison  to  the  ash 
contents  of  terrestrial  biomass  (wood,  grasses,  and  crop  residue) 
which  are  below  10  wt.%.  This  is  explained  by  algae’s  high  potential 
for  mineral  biosorption  and  due  to  their  simple  cellular  structures. 
The  ash  content  in  aquatic  biomass  vary  drastically  from  less  than 
10%  (Gai  et  al.,  2013)  to  more  than  50%  (Zhao  et  al.,  2013),  accord¬ 
ing  to  type  of  alga  and  growth  conditions.  Algae  grown  on  waste- 
water  or  high-mineral  effluents  in  open  ponds/raceways  are 
expected  to  contain  higher  ash  than  algae  produced  in  a  closed- 
system  on  low-minerals  water.  The  difference  in  ash  contents,  12 
percentage  points,  in  this  study  could  be  attributed  to  use  of  the 
alum  flocculant,  A12(S04)314H20,  on  the  swine  manure  slurries 
used  to  grow  Algae  2.  The  high  percentages  of  A1203  and  S03  found 
in  Algae  2  ash  residue,  19.86%  and  10.25%,  compared  to  Algae  1, 
1.93%  and  4.54%,  respectively,  supports  this  hypothesis  (Table  1). 
This  difference  in  the  ash  contents  also  explains  the  relatively 
low  volatile  matter  and  fixed  carbon  contents  observed  in  Algae 
2  when  compared  to  Algae  1 . 


3.3.  Elemental  composition,  and  higher  heating  value 

Percentages  of  the  main  organic  elements:  carbon  (C),  hydrogen 
(H),  nitrogen  (N),  and  oxygen  (O),  in  both  algae  are  shown  in 
Table  1.  Algae  1  sample  showed  higher  content  of  organic  elements 
when  compared  to  Algae  2  sample,  except  for  N.  These  differences 


Table  1 

Proximate,  ultimate,  ash,  calorific  value,  and  pH  analyses  for  the  two  algae  consortia 
studied. 


QOFtg(%)  =  100 


V^Pl(mcalc-m0bs)2 


Nn 


(16) 


Numerical  integration  was  performed  using  a  numerical 
integrator  for  stiff-equations,  (ODE15s)  while  the  constrained  non¬ 
linear  minimization  was  performed  using  the  (fmincon)  function, 
both  of  which  are  part  of  the  MATLAB  software  package  (MATLAB 
R2013b,  The  Mathworks,  Inc.,  Natick,  MA,  USA). 


3.  Results  and  discussions 

3.1.  Species  identification 

The  harvested  microalgae  species  were  inspected  non-quantita- 
tively  using  optical  microscopy.  The  microalgae  consortium  in  each 
harvest,  i.e.,  Algae  1  and  Algae  2,  consisted  mostly  of  filamentous 
microalgae  mixtures  in  addition  to  diatoms.  The  most  common 
filamentous  genus  identified  in  Algae  1  was  Mougeotia,  while  the 
genus  Cladophora  dominated  Algae  2  harvest.  Both  genera  are  com¬ 
mon  fresh-water  microalgae,  which  were  reported  to  be  tolerant  of 
a  wide  range  of  growing  conditions  (Higgins  et  al.,  2008;  Morris, 
1982).  In  a  study  of  algal  growths  in  Lake  Mead  (Nevada-Arizona, 
USA),  Morris  (1982)  also  reported  the  predominance  of  Mougeotia 
species,  in  certain  monitoring  stations  during  the  summer  months 
(June-July)  which  then  changed  to  Cladophora  and  Stigeoclonium 
during  the  months  of  October  and  November. 


Algae  1 

Algae  2 

Proximate  analysis  (wt.%) 

Moisture 

5.14 

5.09 

Volatile  matter 

64.49 

56.30 

Fixed  carbon 

10.41 

6.32 

Ash 

19.96 

32.29 

Ultimate  analysis  (wt.%) 

C 

41.51 

33.79 

H 

5.59 

4.73 

N 

5.40 

6.35 

S 

0.51 

1.57 

O 

27.03 

21.27 

Ash 

19.96 

32.29 

Ash  analysis  (wt.%) 

A1203 

1.93 

19.86 

Si02 

14.57 

16.44 

p2o5 

15.69 

15.27 

I<20 

6.57 

12.09 

so3 

4.54 

10.25 

CaO 

37.15 

9.85 

Na20 

3.08 

4.88 

MgO 

4.46 

4.06 

Fe203 

1.42 

3.50 

MnO 

0.08 

0.17 

Ti02 

0.06 

0.15 

HHV  (kjg  ’) 

16.63 

14.53 

PH 

6.88 

7.69 

Weight-basis  percentage. 
By  difference. 

Higher  heating  value. 
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are  attributable,  as  discussed  earlier,  to  the  use  of  metal  salt  to  pre¬ 
cipitate  the  suspended  solids  in  the  manure  slurry  used  in  Algae  2 
production.  By  contrast,  the  higher  N  content  in  Algae  2  is  due  to 
the  elimination  of  slurry  dilution,  which  was  employed  in  growing 
the  Algae  1  consortia.  N  loading  was  reported  to  influence  the  lipid 
and  protein  in  the  algal  biomass  (Xin  et  al.,  2010).  In  lipid  produc¬ 
tion  from  microalgae,  a  nutrient  stress  is  typically  imposed  by 
reducing  N-loading  to  induce  lipid  production  and  storage  in  algae 
cells.  In  phycoremediation,  by  contrast,  algae  are  grown  on  N-rich 
substrates,  which  increase  the  protein  content  in  the  cells.  From 
their  N  contents,  the  protein  content  in  Algae  1  and  2  can  be 
estimated  using  Jones’  factor  of  6.25  (Jones,  1941)  to  be  33.8% 
and  39.7%,  respectively.  Unlike  terrestrial  biomasses,  which  consist 
primarily  of  cellulose,  hemicellulose,  and  lignin,  microalgae  are 
single  cell  organisms  that  consist  of  carbohydrates,  proteins,  and 
lipids.  This  structural  difference  can  be  observed  in  the  stoichiom¬ 
etric  expressions  for  Algae  1:  CHj  62No.nOo.49  and  Algae  2:  CHj  6g 
N0.i6O0.47,  as  compared  to  that  for  hybrid  poplar  wood  (Jenkins 
et  al.,  1998):  CH1.45N0.01O0.60- 

The  difference  observed  in  the  higher  heating  values  (HHV) 
between  the  studied  algae  harvests,  16.6  kjg-1  in  Algae  1,  and 
14.5  kj  g_1  in  Algae  2,  could  be  attributed  to  the  difference  in  the 
ash  contents.  Calculating  the  dry,  ash-free  HHV  (DAF),  using 
the  ash  contents  and  the  dry-basis  HHVs  in  Table  1,  shows  both 
algae  to  have  similar  energy  contents,  i.e.,  21  kj  g^p. 

3.4.  Thermogravimetric  analysis 

With  the  temperature  increase,  the  samples  underwent  a  series 
of  endothermic  and  exothermic  reactions,  which  involved  varying 
degrees  of  weight-loss  (devolatilization)  and  structural  change  to 
the  sample  matrix.  In  both  algae  consortia,  decomposition  temper¬ 
atures  increased  with  increasing  the  heating  rate  (/)).  Doubling  p 
led  to  an  increase  between  7  and  9  °C  for  the  entire  weight  loss 
curve.  Similarly,  the  rate  of  weight  loss  was  observed  to  double 
with  each  doubling  of  the  heating  rate.  In  Algae  1,  the  pyrolysis 
resulted  in  a  weight  loss  of  69-72%  of  initial  weight  whereas  in 
Algae  2  the  weight  loss  ranged  between  62%  and  67%  over  the 
entire  temperature  range,  30-800  °C.  The  difference  in  ash  con¬ 
tents  between  the  two  algal  consortia  was  reflected  in  difference 
between  overall  weight-loss. 

During  algae  pyrolysis,  the  weight  loss  can  be  divided  into  three 
consecutive  stages:  drying,  active  pyrolysis,  and  passive  pyrolysis. 
Drying,  which  typically  takes  place  below  110  °C,  is  the  first  step  in 
algae  pyrolysis  and  it  involves  the  evaporation  of  both  free  and 
hygroscopic  water  in  the  sample  matrix.  In  this  step  the  initial 
weight  dropped  around  5%  with  the  corresponding  DTG  peak  at 
around  100  °C.  After  drying,  the  samples  underwent  a  brief  induc¬ 
tion  period,  between  110  and  150  °C,  in  which  the  weight  loss  was 
minimal. 

In  pyrolysis,  most  of  the  weight  loss  takes  place,  between  150 
and  550  °C,  during  the  active  pyrolysis  stage.  This  weight  loss  rep¬ 
resents  the  thermal  depolymerization  and  volatilization  of  various 
organic  matter  (volatile  matter)  components.  Algae  1  lost  58%  of 
total  weight  during  the  active  pyrolysis  step  compared  to  only 


50%  weight-loss  in  Algae  2.  Despite  the  similarity  in  peak  decom¬ 
position  temperatures  in  the  two  consortia  (Table  2),  the  maxi¬ 
mum  decomposition  rates  for  Algae  1  samples,  0.44-0.48%  °C_1, 
were  consistently  higher  than  those  observed  for  Algae  2, 
0.33-0.37%  °C-1. 

In  biopolymers  such  as  microalgae,  the  thermal  decomposition 
represents  the  summation  of  decomposition  profiles  for  the  indi¬ 
vidual  components.  Given  the  wide  variability  between  microalgae 
species,  in  terms  of  main  ingredients’  concentration  and  composi¬ 
tion,  the  decomposition  profiles  under  pyrolysis  can  vary  greatly. 
Maddi  et  al.  (2011)  compared  pyrolysis  thermograms  of  Lyngbya 
sp.  and  Cladophora  sp.,  and  attributed  the  dissimilarity  to 
differences  in  protein  compositions  (Maddi  et  al.,  2011).  In  this 
study,  the  seasonal  variability  between  consortia  induced  differ¬ 
ences  in  algal  community  compositions,  which  were  reflected  in 
the  respective  decomposition  profiles.  In  both  consortia,  however, 
the  main  decomposition  peak  took  place  in  the  range  associated 
with  carbohydrates  and  proteins  indicating  a  similarity  in  the  main 
composition.  Furthermore,  the  peak  associated  with  lipids  decom¬ 
position  (above  390  °C)  (Na  et  al.,  2011)  was  only  a  small  shoulder 
in  both  consortia,  indicating  low  lipid  contents. 

During  passive-pyrolysis  (T>550°C),  both  algae  harvests 
showed  minimal  weight-loss,  which  is  characteristic  of  char 
pyrolysis.  A  minor  decomposition  peak  (0.06  %  °C)  which  was 
observed  in  Algae  1  between  650  and  700  °C  was  not  observed  in 
Algae  2.  In  studying  the  pyrolysis  of  Chlorella  vulgaris,  Agrawal 
and  Chakraborty  (2013)  attributed  a  similar  high  temperature 
weight-loss  peak  (at  700  °C)  to  a  char  gasification  step,  which 
was  explained  as  C02  loss  (Agrawal  and  Chakraborty,  2013). 
Maddi  et  al.  (2011)  also  reported  a  high  temperature  peak  during 
the  decomposition  of  Cladophora  sp.  but  could  not  identify  the 
component  associated  with  it.  The  kinetics  of  decomposition 
for  both  algae  consortia  during  the  active  pyrolysis  stage,  150-550  °C, 
will  be  discussed  in  the  following  section. 

3.5.  Pyrolysis  kinetics 

The  temperature-conversion  data  in  each  algae  sample  (Fig.  1) 
was  used  to  determine  the  decomposition  kinetics  during  the 
active  pyrolysis  stage  (150-550  °C).  The  conversion  profile  for  each 
algae  consortium  approximates  a  sigmoidal  function,  which  is 
typically  associated  with  the  autocatalytic  reactions  involved  in 
pyrolysis  (McCoy,  1999). 

The  kinetics  determined  using  isoconversional  methods,  i.e., 
activation  energy  and  pre-exponential  term,  are  detailed  in  Table  3. 
Both  isoconversional  and  integral  methods  showed  strong  correla¬ 
tion  coefficients  (r)  and  high  significance  (p-value  <  0.01).  In  each 
algae  consortium,  the  KAS  and  FWO  activation  energy  values  were 
very  close,  within  2%,  at  each  conversion  degree.  Activation 
energies  determined  using  Friedman  method,  on  the  other  hand, 
were  consistently  higher  than  KAS  and  FWO,  at  times  by  12%,  as 
shown  in  Fig.  2.  This  difference  in  the  activation  energies,  between 
the  differential  and  integral  methods,  could  be  attributed  to  the 
numerical  differentiation  step  necessary  to  determine  the  kinetics 
in  Friedman  method.  Although  Friedman  method  does  not  rely  on 


Table  2 

Pyrolysis  temperature  (°C)  and  weight  loss  (%)  corresponding  to  peak  decomposition  rate  (%  <  '}  in  both  algae  consortia. 


Heating  rate  (°C  min  *) 

Algae  1 

Algae  2 

-(dW/dT)p(%°  C1) 

Tp  (°C) 

WLP  (%) 

-(dW/d  T)P(%“C-1) 

TP  (°C) 

WLP  (%) 

5 

0.48 

298.58 

28.4 

0.36 

302.16 

24.4 

10 

0.47 

308.01 

28.2 

0.37 

309.30 

23.9 

20 

0.44 

315.72 

27.3 

0.33 

317.28 

23.9 

40 

0.45 

328.31 

27.6 

0.33 

328.64 

24.3 
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Fig.  1.  Conversion-temperature  (a-T)  curves  for  (A)  algae  1,  and  (B)  algae  2  under 
different  heating  rates:  fi  =  5, 10,  20,  and  40  °C  min1. 


the  integral  approximation  which  is  crucial  to  the  integral  meth¬ 
ods,  the  numerical  differentiation  used  in  this  method  magnifies 
the  instrument  noise  which  necessitates  data  smoothing  prior  to 


analysis  thus  introducing  a  degree  of  uncertainty.  The  activation 
energies  calculated  using  Kissinger  method  were  213.4  and 
247.8  kj  mor1  for  Algae  1  and  Algae  2,  respectively.  These  values 
were  constant  over  the  entire  range  of  conversion  since  they  repre¬ 
sent  one  data  point,  a-T  pair,  from  each  thermogram.  By  contrast, 
the  activation  energies  computed  using  Friedman,  KAS  and  FWO 
methods  showed  a  strong  dependence  on  the  degree  of  conversion 
(a)  in  each  algae  consortium.  This  dependence,  however,  can  be 
divided  into  three  stages:  a  <  0.20  where  the  activation  energies 
fluctuated  strongly  with  the  conversion,  0.20  ^  a  ^  0.60  where 
the  activation  energies  were  not  strongly  dependent  on  the  con¬ 
version  degree,  and  finally,  0.60  <  oc  where  the  activation  energies 
increased  dramatically  with  the  progress  of  decomposition.  The 
increased  £a  -  a  dependence  in  polymers  was  associated  with  par¬ 
allel  reactions,  each  having  different  activation  energy  (Vyazovkin, 
1996).  Yao  et  al.  (2008)  suggested  that  multistep  decomposition 
mechanisms  are  the  cause  for  increases  £a  at  a  >  0.7  during  the 
pyrolysis  of  natural  fibers  (Yao  et  al.,  2008).  It  is  worth  noting, 
however,  that  the  activation  energy  estimated  using  Kissinger 
method  closely  represents  the  mean  activation  energy  value  using 
Friedman,  KAS  and  FWO  during  constant  activation  energy  zone, 
i.e„  0.20  ^  a  ^  0.60. 

Results  of  pyrolysis  modeling  using  independent  parallel  reac¬ 
tions  (IPR)  are  detailed  in  Table  4  as  well  as  in  Figs.  3  and  4.  The 
1PR  model  closely  captured  the  weight  loss  profiles  observed 
during  active  and  passive  pyrolysis  stages  for  both  algae  harvests 
as  shown  in  Fig.  3.  The  quality  of  fit  (QOFDTG  %)  ranged  between 
2.09%  and  3.31%  indicating  only  minor  deviations  between  calcu¬ 
lated  and  observed  weight  loss  data.  First-order  reaction  model 
appears  to  be  capable  of  describing  microalgae  pyrolysis  reactions 
satisfactorily.  Earlier  studies  used  this  model  to  describe  the 
pyrolysis  kinetics  of  various  biomass  residues  such  as  cardoon 
stems  and  leaves,  and  rice  hulls  (Teng  and  Wei,  1998;  Damartzis 
et  al.,  2011).  The  apparent  activation  energies  for  Algae  1  pyrolysis 
reactions  ranged  from  33.26  to  97.66  kj  mol-1  while  the  apparent 
activation  energies  for  Algae  2  pyrolysis  reactions  varied  between 
33.87  and  97.36  kj  mor1.  Contribution  of  the  modeled  pseudo¬ 
components  to  the  overall  decomposition  rate  is  detailed  in 
Fig.  4.  In  both  consortia,  the  peak  devolatilization  of  each 
pseudo-component  took  place  at  a  different  temperature,  starting 
at  100-120  °C  for  pseudo-component  1(47.00-49.65  kj  mol-1) 
which  represents  the  devolatilization  of  moisture  and  light 
hydrocarbons.  Between  300  and  400  °C,  both  pseudo-components 
2  and  3  showed  overlapping  peak  devolatilization  representing 


Table  3 

The  pyrolysis  activation  energy  (Ea),  and  intercept  term  (in  z)  corresponding  to  different  degrees  of  conversion  (<x)  for  Algae  1  and  Algae  2  using  Friedman,  Kissinger-Akahira- 
Sunose  (KAS),  Flynn-Wall-Ozawa  (FWO),  and  Kissinger  methods. 


(ot) 

Friedman  method 

KAS  method 

FWO  method 

£a  (kj  mor1) 

lnz(s-') 

r 

£a  (kj  mor1) 

Inzts-1) 

r 

£a  (kj  mor1) 

inzts-1) 

r 

Algae  1  0.2 

203.3 

37.319 

1.000 

195.7 

32.461 

1.000 

194.7 

47.098 

1.000 

0.3 

210.1 

37.856 

0.999 

205.1 

33.158 

1.000 

204.0 

47.855 

1.000 

0.4 

213.7 

37.738 

0.999 

208.6 

32.874 

1.000 

207.5 

47.617 

1.000 

0.5 

223.4 

38.709 

0.999 

213.9 

32.984 

1.000 

212.8 

47.770 

1.000 

0.6 

247.8 

42.195 

0.998 

228.8 

34.835 

0.999 

227.2 

49.670 

0.999 

0.7 

296.0 

49.348 

0.995 

267.0 

40.565 

0.996 

263.8 

55.467 

0.996 

0.8 

361.7 

57.759 

0.991 

343.9 

51.735 

0.990 

337.5 

66.742 

0.990 

Kissinger  method 

213.4 

39.815 

0.998 

Algae  2  0.2 

263.1 

50.229 

0.997 

267.2 

48.025 

0.997 

262.7 

62.662 

0.998 

0.3 

265.3 

49.248 

0.997 

264.9 

45.597 

0.998 

260.9 

60.300 

0.998 

0.4 

271.6 

49.269 

0.997 

266.8 

44.580 

0.997 

262.9 

59.334 

0.997 

0.5 

290.4 

51.609 

0.996 

276.7 

45.197 

0.996 

272.6 

60.000 

0.997 

0.6 

345.6 

60.407 

0.994 

313.8 

50.767 

0.994 

308.1 

65.628 

0.995 

0.7 

454.5 

77.656 

0.991 

417.3 

67.603 

0.990 

406.9 

82.544 

0.990 

0.8 

538.3 

87.512 

0.984 

527.9 

82.781 

0.984 

512.6 

97.834 

0.985 

Kissinger  method 

247.8 

46.886 

0.999 

‘  in  A  (s~’). 
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(a) 


(a) 


Fig.  2.  The  activation  energy  (Ea,  kj  mol  *)  corresponding  to  the  different  degrees  of 
conversion  (a)  during  pyrolysis  of:  (A)  Algae  1,  and  (B)  Algae  2  using  the  following 
methods:  Friedman  (FR),  Kissinger  (K),  Kissinger-Akahira-Sunose  (KAS),  and 
Flynn-Wall-Ozawa  (FWO). 

the  bulk  of  the  total  weight  loss  (between  72%  and  92%  of  total 
weight  loss)  (Table  4).  Decomposition  peak  for  reaction  2,  as  shown 
in  Fig.  4,  coincided  with  the  overall  peak  decomposition  observed 
in  both  algae  as  listed  in  Table  2.  The  two  overlapping  peaks: 


Temp.  (°C) 


Temp.  (°C) 

Fig.  3.  Comparison  between  observed  weight  loss  (markers)  and  calculated  weight 
loss  (lines)  using  independent  parallel  reactions  (IPR)  model  for  Algae  1  (A)  and 
Algae  2  (B). 

reactions  2  and  3  in  both  consortia  appear  to  represent  the  decom¬ 
position  of  the  main  algae  ingredients,  i.e.,  protein  and  starch.  The 
decomposition  peaks  of  individual  starch  and  protein  model  sub¬ 
stances:  corn  starch  and  Lysozyme  protein  (Maddi  et  al.,  2011) 
appear  to  resemble  the  overlapped  peaks  in  Fig.  4.  The  fourth 
pseudo-component  showed  peak  decomposition  around  600  °C, 
which  represents  the  passive  pyrolysis  weight  loss. 

The  differences  observed  between  the  apparent  activation 
energies  determined  using  isoconversional  methods,  and  those 


Table  4 

Pyrolysis  kinetics  of  the  two  algae  consortia  as  modeled  using  Independent  parallel  reactions  (IPR)  model. 


Reaction 

Kinetics 

Algae  1 

Algae  2 

5  °C  min-1 

1 0  CC  min  1 

20°Cmin  1 

40  °C  min  ' 

5  °C  min-1 

10  °C  min~' 

20  °C  min  1 

40  °C  min  1 

1 

Wi 

0.11 

0.01 

0.02 

0.02 

0.04 

0.03 

0.05 

0.03 

4,  (s-1) 

3.3E+01 

3.8E+04 

8.6E+04 

5.4E+04 

1.0E+03 

4.3E+03 

1.0E+02 

9.0E+04 

Ea,  (kj  mop1) 

47.00 

49.17 

49.65 

48.57 

47.06 

46.20 

32.89 

49.63 

2 

W2 

0.31 

0.36 

0.39 

0.42 

0.29 

0.30 

0.35 

0.31 

42 (s-1) 

2.1E+06 

3.1E+06 

3.6E+06 

5.3E+06 

2.1E+06 

2.3E+06 

2.7E+05 

4.1E+06 

Ea 2  (kj  mol  ’) 

97.09 

97.66 

97.03 

97.24 

97.36 

96.35 

84.23 

95.86 

3 

w3 

0.41 

0.49 

0.49 

0.50 

0.49 

0.52 

0.43 

0.51 

43  (s-1) 

1.8E+00 

1.1E+00 

1.7E+00 

4.9E+00 

1.8E+00 

1.2E+00 

5.5E+00 

3.6E+00 

Ea3  (kj  mol  ') 

39.59 

33.60 

33.26 

35.77 

39.59 

34.42 

39.42 

33.87 

4 

w4 

0.16 

0.14 

0.11 

0.07 

0.18 

0.15 

0.17 

0.15 

44 (s-1) 

3.5E+01 

1.4E+01 

2.8E+01 

8.1E+01 

3.8E+01 

8.6E+01 

8.5E+01 

8.1E+01 

Ea4  (kj  mop1) 

78.06 

69.37 

69.35 

69.28 

78.03 

80.49 

74.27 

69.10 

QOF* 

DTG  (%) 

3.31 

2.89 

2.80 

2.76 

2.78 

2.09 

2.54 

3.03 

TG  (%) 

1.62 

1.03 

0.61 

0.55 

0.80 

0.75 

0.56 

0.53 

Quality  of  fit. 
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Temp.  (°C)  Temp.  (°C) 

Fig.  4.  Comparison  between  observed  weight  loss  (markers)  and  calculated  weight  loss  (lines)  using  independent  parallel  reactions  (IPR)  model  for  Algae  1  (A)  and  Algae  2  (B) 
with  the  different  heating  rates,  5,  10,  20  and  40  °C  min-1,  represented  by  the  numerals  1,  2,  3  and  4,  respectively. 


determined  using  the  independent  parallel  reactions  (IPR)  model 
are  attributed  in  part  to  the  reaction  model  used.  While  isoconver- 
sional  methods  determine  an  overall  apparent  activation  energy 
without  specifying  a  reaction  model  (model-free),  modeling  the 
pyrolysis  kinetics  necessitate  the  assumption  of  a  reaction  model 
with  multiple  activation  energies  for  the  different  reactions 
involved.  Damartzis  et  al.  (2011)  showed  a  similar  contrast 
between  the  activation  energies  of  cardoon  pyrolysis  determined 
using  isoconversional  methods  and  those  using  IPR  model 


(Damartzis  et  al.,  201 1 ).  They  suggested  that  isoconversional  meth¬ 
ods  are  most  suited  to  qualitative  evaluations  of  the  pyrolysis  pro¬ 
cess,  whereas  models  such  as  IPR  are  more  suitable  quantitatively 
to  study  and  model  the  pyrolysis  process. 

It  is  worth  noting,  however,  that  in  terrestrial  biomasses  such  as 
wood,  and  crop  residues,  the  individual  pseudo-components 
resulting  from  kinetic  modeling,  such  as  IPR,  can  be  easily  associ¬ 
ated  with  original  biomass  ingredients,  i.e.,  cellulose,  hemicellulose 
and  lignin.  However,  given  the  species  heterogeneity  of  the  algal 
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consortia  investigated  here,  not  to  mention  the  complex  structures 
in  individual  algal  species,  it  is  not  possible  to  map  the  modeled 
pseudo-component  to  specific  structural  components  or  interme¬ 
diate  species.  Nonetheless,  coupling  these  pseudo-components  to 
the  evolved  species,  via  spectral  analysis,  can  offer  better  under¬ 
standing  of  the  original  species  undergoing  decomposition. 

4.  Conclusions 

•  The  algae  genus  Mougeotia  was  the  common  genus  in  the 
Algae  1  harvest,  while  the  genus  Cladophora  was  predomi¬ 
nant  in  the  Algae  2  harvest. 

•  In  isoconversional  methods,  the  apparent  activation  ener¬ 
gies  for  pyrolysis  of  Algae  1  were  lower  than  Algae  2 
pyrolysis. 

•  Friedman  method  activation  energy  values  were  within 
12%  of  Kissinger-Akahira-Sunose  (KAS)  and  Flynn-Wall- 
Ozawa  (FWO)  values. 

•  The  activation  energy  dependence  on  conversion  suggests 
complex  reaction  schemes  which  should  not  be  reduced 
to  a  single-step  reaction. 

•  The  pyrolysis  kinetics  of  each  consortium  was  modeled 
using  independent  parallel  reaction  (IPR)  model  as  a  group 
of  four  parallel,  independent,  first-order  reactions. 
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